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1.  INTRODUCTION 


The  use  of  electrostatic  effects  between  droplets  has 
been  postulated  as  a potential  mechanism  for  dissipating 
warm  fog.  Two  separate  effects,  summarized  in  Figure  1, 
must  be  considered  with  respect  to  electrically  enhanced 
coalescence.  One  is  the  re  suit  of  an  applied  electric  field 
which  polarizes  all  droplets.  Tag  (1976)  considered  this 
possibility  and  found  that  extreme  electric  fields  are 
necessary  to  produce  a worthwhile  improvement  in  visibility, 
and  then  only  in  fogs  with  relatively  large  water  droplets. 

The  second  effect  is  the  focus  of  study  in  this  paper. 

The  effect  of  electric  fields  and  charges  on  collision- 
coalescence  has  long  been  recognized.  Lord  Rayleigh  (1879) 
noted  the  effect  of  electricity  on  drops  emitted  from  a water 
jet.  In  1923,  an  article  in  Scientific  American  commented 
on  a proposal  that  charged  sand  particles  be  used  to  clear 
a fog.  Pauthenier  (1950)  suggested  that  charged  water  drops 
be  used.  Carroz  et_  a_l_.  ( 1 972  ) made  an  initial  attempt  at 
determining  visibility  improvement  as  the  result  of  charged 
seedings,  but  considered  only  monodisperse  fog  and  treatment 
drops  in  conjunction  with  a continuous  growth  process  for 
the  charged  drops.  Tag  (1976)  provides  a more  complete 
history  of  both  electrostatic  effects. 

Using  a microphysical  fog  model,  this  study  will  examine 
charged  drop  seeding  as  a fog  dispersal  technique.  The  purpose 
is  to  evaluate  the  potential  of  this  method,  as  well  as  to 
isolate  those  parameters  which  would  most  control  the  degree  of 
visibility  improvement  in  an  experimental  environment. 
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Figure  1.  The  two  separate  aspects  of  electrically 
enhanced  coalescence. 
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2.  PRELIMINARY  CONSIDERATIONS 


Before  proceeding,  a discussion  of  drop  charge  and 
associated  collision  efficiencies  is  warranted.  In  the 
situation  simulated  previously  (Tag,  1976),  a uniform  elec- 
tric field  was  applied  continuously  to  the  entire  fog, 
resulting  in  collision  efficiencies  which  remained  unchanged 
for  any  drop  size-pair  combination.  The  numerical  simulation 
becomes  more  complicated  when  individual  drop  charges  must 
be  conserved  with  every  collision  and  also  result  in  collision 
efficiencies  that  change  for  every  size  combination  with 
eve ry  time  step.  Difficulties  arise  in  finding  an  efficient 
method  of  determining  appropriate  collision  efficiencies. 

a . Drop  Charging 

It  is  necessary  to  understand  the  charging  limita- 
tions appropriate  for  water  drops.  Lord  Rayleigh  (1882)  was 
the  first  to  derive  the  limit  to  which  a water  drop  can  be 
charged  before  rupturing: 

Q = Vl6  7 a r3  eSU  0) 

where  a is  the  surface  tension  and  r the  radius  of  the  drop. 
Equation  (1)  is  commonly  referred  to  as  the  "Rayleigh  Limit" 
and  defines  that  charge  at  which  electrical  forces  which 
tend  to  rupture  the  drop  exactly  balance  the  surface  tension 
holding  the  drop  together.  Empirical  verification  of  this 
limit  has  been  attempted  by  several  workers.  Using  charged 
oil  drops  with  radii  0.1  to  10  microns,  Hendricks  (1962) 
found  very  close  agreement  to  the  Rayleigh  theoretical  curve. 
More  recent  verification  was  made  by  Schweizer  and  Hanson 
(1971).  Using  drops  made  of  n-octanol  with  a size  range  7.5 
to  20  microns,  the  measured  charge  to  mass  ratios  of  the 
drops  at  breakup  were  within  a standard  deviation  of  M of 
the  Rayleigh  Limit. 
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A second  point  that  must  be  considered  is  the  length  of 
time  a drop  can  hold  its  charge.  If  a discharge  should  be 
instantaneous,  no  advantage  would  be  gained  from  an  initially 
large  charge.  Vonnegut  (Moore,  1973),  however,  states  that 
objects  supported  in  the  atmosphere  will  lose  charge  at  a 
rate  corresponding  to  a half-life  of  10  to  20  minutes. 

Schweizer  and  Hanson,  in  the  experiments  described  above, 
found  that  their  drops,  which  were  often  charged  completely 
to  the  Rayleigh  Limit,  occasionally  lost  some  charge,  but 
that  this  loss  was  never  very  large.  In  this  paoer,  it  is 
assumed  that  such  leakage  would  likewise  be  small  in  an  air- 
flow, as  found  in  the  wake  of  an  aircraft.  Because  the 
actual  charges  that  are  simulated  are  only  a fraction  of  the 
Rayleigh  Limit  maximum  (see  below),  and  also  because  of  the 
rapidity  with  which  these  drops  should  grow,  the  assumption 
of  negligible  charge  loss  is  probably  acceptable. 

Finally,  assuming  that  a large  charge  is  desirable  and 
will  remain  on  the  drop,  what  are  the  operational  limitations 
in  achieving  that  charge?  Carroz  and  Keller  (1976)  conducted 
laboratory  experiments  to  isolate  parameters  important  for 
achieving  maximum  drop  charge.  Using  induction  charging, 
they  found  that  the  degree  of  charging  is  a function  of  nozzle 
type,  liquid  flow  rate,  airflow  speed  past  the  nozzle,  voltage, 
and  several  other  factors.  In  general,  they  found  that,  using 
cone-type  nozzles,  drops  could  be  charged  to  a maximum  of  25% 
of  the  Rayleigh  Limit,  with  an  average  value  closer  to  15%. 
These  figures  provide  realistic  limitations  to  the  drop 
charges  that  are  used  in  the  numerical  experiments  conducted 
in  this  study. 

b . Collision  Efficiencies 

In  order  to  ascertain  the  effect  of  charged  drops 
in  a microphysical  model,  one  must  have  precise  values  of 
collision  efficiency  for  the  drop  interactions.  As  outlined 
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in  Tag  (1976),  various  theoretical  and  experimental  studies 
have  defined  efficiencies  for  various  alignments  of  electric 
fields  and  drop  charge  combinations.  Because  of  non-1 i neari ty 
of  the  results,  however,  analytic  solutions  and  numerical 
approximations  have  been  few.  Using  a trajectory  model  based 
on  work  by  Davis  (1964),  and  Davis  and  Sartor  (1967),  Paluch 
(1970)  compared  computed  efficiencies  to  analytically  derived 
values  for  certain  simplified  cases.  Unfortunately , only 
oppositely  charged  drops  were  considered.  For  the  situation 
here,  a charged-neutral  interaction  is  the  one  of  interest. 

Although  not  considered  by  Paluch,  the  trajectory  model 
noted  above  is  capable  of  simulating  charged-neutral  combina- 
tions. Dr.  Doyne  Sartor  of  NCAR  offered  his  computer  program 
for  computing  these  efficiencies.  It  was  intended  that  this 
model  be  used  in  conjunction  with  the  microphysical  fog 
model.  It  soon  became  clear,  however,  that  the  models  could 
not  be  used  together.  On  a CDC  6500,  one  collision  efficiency 
calculation  takes  on  the  order  of  20  seconds.  Since  several 
hundred  efficiency  calculations  would  be  necessary  per  time 
step,  it  was  realized  that  a simplifying  assumption  would  be 
necessary. 

An  earlier  paper  by  Cochet  (1952)  provided  this  simpli- 
fication. Cochet  was  interested  in  the  possibility  of  seeding 
neutral  clouds  with  charged  drops  with  the  idea  of  initiating 
precipitation.  In  the  process  of  calculating  the  growth  of 
a charged  drop,  he  derived  an  expression  for  collision 
efficiencies  for  charged-neutral  combinations.  As  Paluch  did 
in  her  work,  Cochet  made  certain  approximations  in  his 
derivation  from  the  equation  of  motion,  the  primary  one  being 
that  the  inertia  of  the  charged  collector  modifies  negligibly 
the  flow  field  around  itself.  Such  an  assumption  is  quite 
valid  for  small,  highly-charged  collectors,  but  becomes 
increasingly  suspect  as  the  collector  increases  in  size  or 
as  the  charge/mass  ratio  decreases.  Cochet  found  this  approxi- 
mation to  be 
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where  Y£  is  the  linear  collision  efficiency,  x = (e-l)/(e+2) 
and  e the  dielectric  constant,  p the  density  of  water,  q 
charge  (esu),  and  r.  and  r are  the  radii  of  the  collector 
and  collected  drops,  respectively.  Notice  that  collision 
efficiency  is  not  a function  of  the  charge  sign.  In  this 
study,  all  drop  charges  will  be  assumed  to  be  of  the  same 
sign,  and  positive. 

Cochet  recognized  Eq . (2)  as  an  approximation  and 

attempted  to  determine  its  inaccuracy.  Figure  2,  taken 
directly  from  Cochet's  1952  paper,  shows  a comparison  between 
Eq.  (2)  (dot-dashed  lines)  and  actual  values  (solid  lines) 
for  a collector  drop  charge  of  4 X 10  ^ esu  (1  coulomb  = 

3 X 10^  esu).  This  charge  corresponds,  approximately,  to  20 
and  2%  of  the  Rayleigh  Limit  for  10  and  50  microns  radius 
respectively.  The  term  "actual"  refers  to  efficiencies 
obtained  without  the  approximations  inherent  in  the  deriva- 
tion of  Eq . (2).  Cochet  determined  "actual"  collision 
efficiencies  by  performing  trajectory  calculations  graphic- 
ally, instead  of  numerically  as  done  in  the  NCAR  model.  To 
ascertain  the  accuracy  of  his  graphical  technique,  several 
points  on  Cochet's  curves  were  recalculated  using  the  NCAR 
model.  Remarkably,  nearly  all  calculated  values  lay  within 
the  readable  error  in  Figure  2,  due  to  the  care  with  which 
Cochet  drew  his  trajectories  or,  alternately,  to  the  credit 
of  the  NCAR  model.  For  a given  charge,  therefore,  Figure  2 
shows  that  the  approximation  is  generally  an  over-estimator 
as  the  collector  increases  in  size  or,  as  noted  above, 
inertia  becomes  more  and  more  a factor. 
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Figure  2.  Charged-neutral  collision 

efficiencies  (y  --  Cochet's  symbol  for 
this  variable) , actual  (solid)  and  from 
Eq.  2 (dot-dashed) , for  collector 
diameters  D and  various  collected 
diameters,  d,  for  a collector  charge  of 
4.0E-4  esu  (taken  from  Cochet  (1952)). 
The  dashed  lines  represent  hydrodynamic 
efficiencies  (after  Langmuir) . 
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Although  obviously  inaccurate  under  certain  conditions, 
Cochet's  equation  is  used  in  this  study  for  the  following 
reasons.  Firstly,  no  better  approximation  could  be  found. 
Secondly,  the  equation  is  accurate  under  those  conditions 
which  are  being,  at  least  initially,  simulated  (i.e.,  highly 
charged,  relatively  small  droplets).  Because  Figure  2 was 
drawn  for  drops  all  having  the  same  absolute  charge,  more 
error  appears  at  the  larger  sizes  than  would  have  if  all 
drops  had  a certain  percentage  of  the  Rayleigh  Limit  charge. 
Additionally,  the  approximation  is  always  most  accurate  when 
collision  efficiencies  are  largest,  and  thus  when  they  have 
their  most  impact.  Inaccuracies  can  be  better  tolerated  when 
their  effect  is  less.  Finally,  recognizing  that  inaccuracies 
will  almost  always  appear  as  over-estimates,  one  has  an 
understanding  of  the  bias. 
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3.  MODEL  DESCRIPTION 


The  basic  model  is  the  same  as  that  described  in  Tag 
(1976).  At  each  of  ten  model  levels  a fog  droplet  spectrum  is 
specified  and  evolves  as  the  result  of  collision-coalescence, 
sedimentation,  and  turbulent  mixing.  Visibility  is  computed 
from  the  drop  spectra.  The  change  required  in  the  model  for 
this  study  is  the  addition  of  a second  "treatment"  spectrum, 
subject  to  the  same  processes  as  above,  but  which  also  inter- 
acts with  the  foq  spectra.  Because  seedinq  drops  will  carry 
a charqe  and  fog  droplets  are  assumed  neutral,  a separation 
of  the  spectra  is  necessary.  In  addition,  the  total  charqe 
assigned  to  the  treatment  categories  must  be  conserved  and 
fol 1 owed . 

Letting  Np  and  N j (as  a function  of  radius  r and  height  z) 
be  the  average  number  density  of  fog  and  treatment  (charqed) 
drops,  respectively,  the  two  spectra  evolve  according  to  the 
following: 
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where  collision-coalescence  is  represented  by  the  following 
interaction  combinations: 


FF  = Foq  - Foq 

TT  = Treatment  - Treatment 

FT  = Fog  - Treatment 
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Because  two  separate  spectra  are  represented  simultaneously, 
these  three  combinations  are  necessary.  However,  after  a 
few  preliminary  computer  experiments,  it  was  decided  to 
eliminate  FF  interactions.  Since  it  had  been  intended  that 
empirical  spectra  from  mature  fons  be  used,  little  advantaqe 
v/ as  seen  in  allowing  the  fog  spectrum  to  evolve  further. 
Consequently,  FF  interactions  were  set  to  zero. 

Notice  that  there  are  no  fog-treatment  interactions  which 
result  in  a "gain"  for  the  fog  spectrum.  This  convention  was 
chosen  so  as  to  result  always  in  the  dilution  of  the  treatment 
drop  and  its  subsequent  placement  back  into  the  treatment 
spectrum.  This  segregation  is  mandatory  since  the  electro- 
static efficiencies  provided  by  Eq.  (2)  are  valid  only  for 
charged-neutral  combinations.  For  this  reason,  the  resulting 
drop  of  a collision  involving  a treatment  drop  is  always 
regarded  as  a treatment  drop.  Alteration  of  the  charge/mass 
ratio  of  treatment  drops,  as  a result  of  collisions  with 
neutral  and  other  charged  drops,  is  essential  in  determining 
appropriate  collision  efficiencies.  In  addition,  besides 
forcing  segregation  of  the  two  spectra,  this  convention  permits 
movement  of  water  from  one  spectrum  to  the  other  (a  measure  of 
the  rapidity  with  which  treatment  is  working)  to  be  followed 
quite  easily. 

The  above  terms  involving  collision-coalescence  can  be 
summarized  by  the  following: 

(5) 

i - 1 M 

-NT( i ) E PT(i »n)B(i ,n)  - E NT(n)PT(n,i)  = TT ( i )Loss  jAt 
T ri  = 1 n = i +1  1 


i -2 

N T ( i - 1 ) E M1"1  » n ) B ( i - 1 ,n) 
1 n=  1 1 


TT  ( i ) 


Gain 


(6) 
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^ 1 ^ Loss  FAt 


when  A=  F , B = T 


when  A=T,B=F 


(7) 


i-2 


i -2 


(8) 


As  in  the  first  part  of  this  study,  B is  the  drop  partitioning 
factor,  i and  n refer  to  category  numbers  and  M is  the  total 
number  of  categories.  P(i,n)  is  defined  as  the  probability 
that  an  ith  will  collect  an  nth  category  drop  in  time  At: 


where  V(i,n),  as  defined  by  Berry  (1967),  is  the  collection 
kernel : 


Yc  is  the  linear  collision  efficiency  as  used  by  Shafrir  and 
Neiburger  (1963)  and  E is  the  coalescence  (as  opposed  to 
collision)  efficiency.  For  reasons  similar  to  those  outlined 
in  Tag  (1976),  the  coalescence  efficiency  will  be  assumed  to 
be  unity.  The  mathematical  approximations  of  Wobus  e^t  al  . 
(1971)  are  used  for  the  drop  terminal  velocities,  w(r). 


P ( i , n ) = V ( i , n ) N(rn)  At 


(9) 


V ( i , n ) = tt  r . ^ 


[Vc(ri  ’rn'q1 }]2  E[w(ri)  ' w(rn)]  (,0) 
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Assuming  that  each  treatment  drop  is  initially  assigned 
an  electrostatic  charge,  the  movement  of  this  charge  as  the 
result  of  collision-coalescence  and  vertical  transfer  must  be 
followed.  Several  methods  of  charge  allocation  and  conserva- 
tion were  considered.  It  was  decided  that  a continuing 
assignment  of  a single  charge  per  drop-category  would  be 
adequate,  rather  than  maintaining  an  array  of  charges  per 
category  as  collection  proceeds.  The  justification  for  this 
decision  stems  from  Eq.  (2).  An  analysis  of  Cochet's  equation 
shows  that  variation  of  collector  drop  radius  is  much  more 
important  than  variation  of  charge.  A doubling  of  radius 
results  in  an  order  of  magnitude  change  in  collision  effi- 
ciency; a similar  doubling  of  charge  results  only  in  a 
doubling  of  efficiency.  For  this  reason,  assignment  of  an 
average  charge  per  drop  category  is  a reasonable  assumption. 
Consequently,  for  collection,  the  procedure  involves  only 
simple  averaging  of  charge  for  those  drops  entering,  as  well 
as  those  remaining  in,  a given  size-category: 


q-j-  ( i ^ OLD  ( ^ j ( 1 ) + A N y ( i ) ) + q jj ( i ) TT ( i ) 


GAIN  T + qFT(i  )FT^~  ^ GA I N T 


rT(i ^GAIN  T 


where , 


(11) 


ANt(1 ) 


tt^loss  t + ft^)loss  t and’ 


i-2 

qTT(i)  = N T ( i - 1 ) £ PT ( i - 1 ,n)B(i-l  ,n) 

1 1 n=  1 


TT^  ^ GAIN  T At 
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and , 


i -2 

qFT(i)  = Np(i-l)  2]  PT ( i - 1 ,n)B(i-l ,n)  qT(n) 


FT^  ^ GAI N TAt’ 


(13) 


recalling  that  q p ( 1 - 1 ) and  qp(n),  by  definition,  are  zero. 

A similar  procedure  is  used  for  conservation  of  charge  as  the 
result  of  movement  between  model  levels.  Although  Eqs.  (11)- 
(13)  conserve  total  drop  charge,  they  more  importantly  permit 
the  charge-to-mass  ratio  to  change  as  the  result  of  collection. 
It  is  this  ratio  which  determines  the  rate  of  collection  of 
fog  droplets. 

In  assigning  appropriate  collision  efficiencies  for  the 
several  collection  combinations,  it  will  be  assumed  that  only 
treatment-collecting-fog  interactions  have  collision  effi- 
ciencies other  than  the  normal  hydrodynamic  ones.  Cochet's 
equation  is  valid  only  for  charged-neutral  collections. 
Obviously,  the  assumption  that  hydrodynamic  efficiencies  apply 
for  charged  collecting  charged,  and  neutral  collecting  charged, 
droplets  is  incorrect.  Because  no  simple  approximation  is 
available  for  these  cases,  hydrodynamic  efficiencies  are 
arbitrarily  used.  Two  considerations,  however,  minimize  the 
effect  of  this  assumption.  First,  although  all  treatment 
drops  of  different  size  have  a different  charge  (a  percentage 
of  the  Rayleigh  Limit),  it  is  assumed  that  all  charges  are 
of  the  same  sign.  Consequently,  collisions  should  not  be 
enhanced  for  charged-charged  interactions.  Second,  there 
will  be  few  FT  collisions  in  which  the  collector  is  a fog 
drop  since  the  charged  drops  must  generally  be  larger  than  the 
fog  droplets  for  the  simulated  mechanism  to  be  effective. 
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4.  AN  ANALYSIS  OF  THE  COLLECTION  KERNEL 


Prior  to  conducting  the  numerical  simulations,  a look  at 
the  collection  kernels  (V)  for  charged  and  uncharged  collectors 
is  warranted.  The  kernel  provides  a much  better  idea  of 
sweepout  effect  than  the  collision  efficiency  alone  since  the 
fall  velocity  difference  is  included  in  V.  Figure  3 gives 
a comparison  between  the  neutral  (Figure  3a)  and  two  charged 
collector  cases.  The  numerical  fits  by  Berry  (1967)  of  the 
Shafri r-Nei burger  collision  efficiencies  are  used  for  the 
neutral  collectors,  and  Cochet's  equation  provides  electro- 
static efficiencies  for  tne  15  and  30%  Rayleigh  Limit  cases. 
These  percentages  provide  reasonable  outlines  to  the  experi- 
mental limits  cited  earlier.  Note  that  the  absolute  collector 
charge  changes  with  radius  since  the  Rayleigh  Limit  increases 
with  radius. 

A quick  comparison  of  the  three  cases  reveals  the  strong 
effect  of  charge  on  the  sweepout  rate.  Increases  from  one 
to  two  orders  of  magnitude  are  evident.  A closer  examination 
shows,  however,  that  such  increases  are  sharply  curtailed 
as  the  collector  size  increases  beyond  25  microns!  As  a 
result,  one  would  conclude  that  seeding  drops  should  be 
relatively  small  in  size.  This  conclusion  is  enhanced  when 
one  considers  that  the  volume  of  liquid  required  increases 
as  the  radius  cubed . 

A comparison  of  Figures  3b  or  c to  3a  shows  that,  in 
addition  to  increasing  in  absolute  magnitude,  the  curve  for 
any  size  collector  is  shifted  to  the  right  --  implying  better 
collection  of  nearly  equal-sized  drops.  This  shift  results 
from  the  rapid  increase  of  collision  efficiency  as  the 
collected  drop  size  approaches  the  collector  radius. 
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Figure  3.  Collection  kernels  for 
several  neutral  (a)  and  charged 
collectors  (b  and  c for  15  and 
30%  of  the  Rayleigh  Limit, 
respectively) . Dashed  lines 
represent  collected  droplet  radii 
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One  additional  point  is  revealed  by  the  dashed  lines  in 
Figures  3b  and  c.  These  lines  follow  two  sizes  of  collected 
droplets,  5 and  10  microns  radius  --  two  sizes  typical  of 
fog  droplet  spectra.  Both  for  3b  and  c,  the  larger  fog  radius 
has  the  better  chance  of  being  collected.  This  observation 
again  implies,  as  was  also  the  case  for  an  applied  electric 
field,  that  results  will  be  closely  tied  to  the  type  of  fog 
being  treated. 
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5.  MODEL  INPUT 


The  two  fog  types  used  in  Tag  (1976)  are  reused  here  and  are 
shown  in  Figure  4.  The  measurements  responsible  for  Fog  A 
were  taken  in  the  Panama  Canal  Zone  and  are  typical  of  fogs 
forming  near  a maritime  influence.  This  spectrum  has  a mean 
volume  radius  of  10.8  microns,  a liquid  water  content  (LWC) 
of  0.39  gm  m”3,  and  produces  a calculated  visibility  of  143  m. 

Fog  B is  based  on  data  taken  in  the  valley  fogs  of  California 
in  which  LWC's  are  typically  below  0.1  gm  m"3,  maximum  drop 
radii  are  approximately  7-8  microns,  and  mean  sizes  average 
around  3 microns  (Reinking,  1975).  Fog  B results  in  a cal- 
culated mean  radius  of  2.9  microns  and  a visibility  of  77  m. 
Spectra  from  actual  data  from  valley  fogs  yielded  spectral 
shapes  very  similar  to  that  of  Fog  A;  consequently,  the  same 
spectral  shape  is  used  for  both  fogs.  These  two  cases  outline 
realistic  limits  of  several  important  fog  parameters  and  should 

allow  assessment  of  these  limits  with  regard  to  seeding.  As 

2 - 1 

in  Tag  (1976),  a constant  exchange  coefficient  of  2 m sec 
is  again  used. 

Probably  the  most  important  input  data  are  the  seeding 
spectrum  and  its  associated  charge.  Carroz  and  Keller  (1975), 
after  assessing  the  advantages  of  various  nozzle  types  for 
an  operational  seeding  program,  chose  30°  cone  angle  nozzles 
with  4.2  gallon  per  hour  (GPH)  flow  rates  (at  300  psi). 

Limited  data  were  available  for  this  nozzle.  However,  a size 
distribution  was  available  for  a 1.0  GPH  nozzle  and  the  spray 
patterns  are  considered  similar  for  the  two.  The  tabular 
portion  of  Figure  5 represents  empirical  data  for  this  nozzle 
(Tate  and  Olson,  1962);  the  accompanying  plot  shows  the  mass 
distribution  for  a 1 gm  m"3  total  water  mass.  A charge  is 
assigned  to  each  size  category  in  the  model  appropriate  to 
this  seeding  distribution.  Using  15%  of  the  Rayleigh  Limit 
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Figure  4.  Water  mass  distribution  for  Fogs  A and  B. 
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Figure  5.  Water  mass  distribution  for  treatment  spray 

spectrum.  Tabular  inset  is  empirical  data  from  which  the 
spectrum  was  fashioned  (after  Tate  and  Olson  (1962)). 
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as  an  average  charge,  Table  1 lists  this  charge  for  each  of 
the  model's  size  categories  applicable  to  the  distribution 
shown  in  Figure  5. 


Table  1 


Radius  (ym) 

Charge  (esu) 

Radius  (ym) 

Charge  (esu) 

4.44 

8.49  E-5 

19.7 

7.93  E-4 

5.23 

1 .08  E-4 

23.2 

1.01  E-3 

6.18 

1.39  E-4 

27.4 

1 . 30  E-3 

7.29 

1 . 79  E-4 

32.3 

1.67  E-3 

8.60 

2.29  E-4 

38.1 

2.13  E-3 

10.1 

2.91  E-4 

45.0 

2.74  E-3 

12.0 

3.77  E-4 

53.1 

3.51  E-3 

14.1 

4.80  E-4 

62.7 

4.50  E-3 

16.7 

6.19  E-4 

73.9 

5.76  E-3 

As  in  the  electric  field  simulations,  a Az  of  15  m is 
used,  which  results  in  a fog  of  150  m depth.  This  depth  is 
realistic  for  a coastal  fog  and  not  uncommon  for  valley  fogs. 
The  15  m grid  spacing  was  chosen  for  two  additional  reasons. 
First,  the  150  m provides  a realistic  height  above  ground 
from  which  aircraft  operations  could  be  safely  carried  out. 
Second,  the  aircraft  from  which  seeding  operations  have,  in 
practice,  been  planned  is  the  B - 2 6 (Reinking,  1975  ).  An 
empirical  estimate  of  the  vertical  diffusion  of  a seeding 
plume  injected  into  this  aircraft's  wake  is  25-30  m (see 
below).  Insertion  at  two  model  levels  provides  a reasonable 
simulation  of  the  vertical  depth  through  which  material  would 
be  distributed. 

It  is  necessary  to  consider  both  the  type  of  aircraft 
and  its  speed  in  determining  the  concentration  of  water  that 
is  to  be  distributed  into  the  aircraft  wake.  Smith  and 
MacCready  (1963)  state  that  the  flow  field  of  an  aircraft 
wake  is  contained  in  an  oval  area  A = 1.69  b2,  where  b is 
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the  aircraft  wingspan.  Assuming  an  even  distribution  within 
this  wake.  Figure  6 relates  seeding  concentrations  to  seeding 
rates  for  several  speeds  of  the  B-26  aircraft  (b  = 21.79  m). 
While  there  are  sure  to  be  wide  variations  of  concentration 
within  the  plume,  the  use  of  this  diagram  does  provide  a 
first  order  estimate  of  liquid  concentration.  Reinking 
(1975),  in  describing  an  airborne  spray  system  for  charged 
drops,  stated  that  the  system  was  capable  of  spraying  at  a 
rate  of  55  gal  min"1.  Depending  on  the  speed  of  the  aircraft, 
Figure  6 suggests  that  this  rate  results  in  concentrations 

. . _ 3 

ranging  from  0.05  to  0.10  gm  m ; such  estimates  will  provide 
a basis  for  model  seeding  concentrations. 
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SEEDING  CONCENTRATION  (jm  m'3 


Figure  6 . 
seeding 


Seeding  concentrations  resulting 
rates  of  the  B-26  aircraft. 


from  various 
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6.  NUMERICAL  RESULTS 


a . Experiment  Set  I - Variation  of  Fog  Type 

Using  the  fog  and  seeding  spectra  described  in  section 
5,  six  numerical  experiments  were  completed.  Two  of  these 
were  designated  as  controls  and  represented  water  seedings 
with  no  associated  charge.  All  experiments  were  run  for 
15  minutes  fog  time  with  a time  step  of  2.5  sec.  It  was 
felt  that  this  period  of  time  provided  a reasonable  goal  for 
operational  usefulness.  In  addition,  fog  regeneration  and 
mixing  into  the  seeded  area  (effects  for  which  the  model  has 
no  capability)  would  greatly  mask  results  beyond  this  time 
period. 

Using  the  nomogram  shown  in  Figure  6,  a 55  gal  min-^ 

seeding  rate  can  produce  a maximum  treatment  concentration 

_ 3 

of  approximately  0.1  gm  m . The  numerical  seedings  of  Fogs 
A and  B were  made  using  this  concentration  and  the  spectrum 
described  in  Figure  5.  Assuming  an  aerial  release  at  the  top 
of  the  fog,  material  is  inserted  at  levels  9 and  10  -- 
corresponding  to  the  depth  suggested  previously.  To  assess 

the  potential  of  a larger  seeding  rate,  the  remaining  two 

_ 3 

cases  were  run  with  0.2  gm  m concentrations. 

Figure  7 gives  visibility  improvement  as  a function  of 

height  for  the  above  seeding  variations  in  Fog  A (initial 

visibility  = 143  m).  Figure  7a  represents  seeding  with 

el ectri cal ly-neutral  water  drops,  using  the  larger  seeding 

_ 3 

concentration  of  0.2  gm  m . Although  a small  change  does 
occur,  most  of  the  improvement  at  fog  top  is  only  the  result 
of  preferential  fallout  of  the  larger  fog  particles.  The 
neutral  seeding  drops  are  relatively  inefficient  collectors. 
Figure  7b  represents  seeding  drops  charged  to  15%  of  the 
Rayleigh  Limit,  as  given  in  Table  1.  Figure  7c  represents 
the  same  seeding  spectrum  but  with  twice  the  seeding  amount. 
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Figure  7 . Visibility  as  a function  of  time  (min)  and  height  for  three  seeding 
variations  in  Fog  A. 


The  effect  of  charge  on  visibility  improvement  is 
evident.  Maximum  impact  on  visibility  occurs  at  the  top  of 
the  fog  where  the  charge/mass  ratio  of  the  seeding  drops  is 
largest.  As  collection  of  neutral  fog  droplets  proceeds,  the 
charge/mass  ratio  of  these  collectors  decreases.  Counter- 
acting this  effect,  however,  is  the  fact  that  the  collectors 
are  increasing  in  absolute  size.  Consequently,  the  effect  of 
charge  density  decrease  is  offset  somewhat  as  the  seeding  drops 
fal 1 through  the  fog . 

As  shown  in  Figure  7c,  a larger  seeding  concentration 
re  suits  (at  least  for  this  case)  in  more  visibility  improve- 
ment. The  transfer  of  water  from  the  fog  to  the  seeding 
spectrum  for  this  experiment  is  demonstrated  in  Figure  8.  At 
alternate  levels  of  the  model,  LWC  for  fog  and  treatment  is 
plotted  as  a function  of  time.  Seeding  material  has  an  initial 
maximum  at  135  and  150  m (level  10  not  shown).  As  water  moves 
from  one  spectrum  to  the  other  as  the  result  of  collection, 
this  maximum  moves  downward  at  a speed  corresponding  to  the 
fall  velocity  of  the  average  drop  sizes.  More  than  80%  of  the 
fog  LWC  is  removed  at  all  levels.  A look  at  the  spectra, 
however,  indicates  that  water  is  preferentially  removed  at  the 
larger  sizes.  Whereas  over  90%  of  20  pm  droplets  are  removed 
at  15  m,  only  45%  of  3 pm  fog  particles  are  swept  out.  This 
result  is  consistent  with  the  collection  kernel  in  Figure  3b. 

This  bias,  unfortunately,  favors  the  wrong  end  of  the  spectrum. 

_ 3 

An  LWC  of  0.1  gm  m which  at  20  pm  radius  will  reduce 
visibility  to  less  than  1000  m,  will  restrict  visibility  at 
3 pm  to  approximately  50  m (Johnson,  1972).  For  this  reason, 
the  80%  LWC  reduction  is  not  as  impressive  as  it  might  be. 

Figure  9 represents  indentical  seedings  of  Fog  B,  the 
narrow  spectral,  small  droplet  fog  typical  of  inland 
California.  The  calculated  initial  visibility  is  77  m. 
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LIQUID  WATER  CONTENT  |gm  nr1) 


Figure  8.  For  alternate  model  levels,  the  change  in 
fog  and  treatment  LWC  for  experiment  depicted  in 
Figure  7c  (Fog  A with  charged  seeding  of  0.2  gm  iri  ) . 
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Figure  9.  Visibility  as  a function  of  time  (min)  and  height  for  three  seeding 
variations  in  Fog  B. 


Examination  of  the  visibility  response  shows  a considerable 
decrease  in  effect  --  noting  that  the  scale  of  the  abscissa 
has  been  changed.  Defining  a visibility  improvement  factor 
(VIF)  as  the  visibility  after  15  min  divided  by  the  initial 
visibility,  the  VIF  decreases  by  more  than  a factor  of  two 
in  going  from  Fog  A to  Fog  B.  Such  a decrease  is  consistent 
with  collection  kernel  ratios  (noted  earlier)  using  average 
radii  for  the  two  fog  spectra  (11  vs  3 ym)  . Figure  10, 
showing  the  time  change  of  fog  and  treatment  LWC's,  indicates 
the  smaller  percentage  (50%  at  15  m)  of  fog  water  removal.  A 
decrease  in  the  downward  speed  of  the  treatment  water  maximum 
indicates  further  the  reduction  of  treatment  drop  growth.  As 
was  also  the  case  for  Fog  A,  removal  of  fog  droplets  increases 
with  increasing  fog  droplet  size,  but  percentages  have 
decreased  considerably. 

b . Experiment  Set  II  - Variation  of  Seeding  Drop  Size 

Experiment  Set  I established  that  the  type  of  fog  as  well 
as  the  quantity  of  seeding  material  will  affect  the  degree  of 
visibility  improvement.  As  is  evident  from  Figure  3,  the  size 
of  the  seeding  drops  is  equally,  if  not  more,  important  than 
the  size  of  the  collected  fog  particles.  The  larger  the 
collector  drop  the  better.  Unfortunately,  an  airborne  spray 
system  imposes  limits  on  both  the  seeding  rate  and  the  quan- 
tity of  available  water.  However,  because  Figure  3 is  plotted 
"per  drop"  rather  than  normalized  per  unit  of  collector  water 
volume,  it  fails  to  show  that  a halving  in  collector  radius 
increases  drop  number  density  by  a factor  of  eight  (for  the 
same  LWC).  Consequently,  although  the  sweepout  rate  of  a 
single  smaller  drop  is  less,  the  combined  effect  of  several 
drops  of  that  size  can  equal  or  exceed  one  larger  drop.  One 
could  thus  hypothesize  that  more  and  more  smaller  particles 
would  be  ever  better.  However,  because  the  final  terminal 
velocity  of  the  collector  drop  should  be  large  enough  to 
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Figure  10.  For  alternate  model  levels , the  change  in 
fog  and  treatment  LWC  for  experiment  depicted  in 
Figure  9c  (Fog  B with  charged  seeding  of  0.2  gm  m 3) . 
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allow  it  to  fall  from  the  fog,  there  must  be  a limit  to  which 
drop  size  can  be  reduced. 

Several  numerical  experiments  confirmed  that  a treatment 

_ 3 

LWC  of  0.1  and  0.2  gm  m at  larger  mean  sizes  (than  that  of 
the  spectrum  in  Figure  5)  produced  less  visibility  improvement. 
To  determine  if  smaller  sizes  would  be  preferable,  the  cases 
described  in  Set  I were  repeated  with  a second  seeding 
spectrum.  All  drop  sizes  shown  in  Figure  5 were  arbitrarily 

divided  by  two  (hereafter  referred  to  as  Seeding  Spectrum  II). 

_ 3 

LWC's  were  normalized  to  0.1  and  0.2  gm  m as  before,  and  a 
charge  of  15%  of  the  Rayleigh  Limit  was  assigned  to  each 
category.  The  22.2  ym  mean  radius  of  Seeding  Spectrum  I was 
thus  changed  to  11.1  ym  in  the  new  drop  distribution. 

Figures  11  and  12  show  visibility  changes  for  Fogs  A and 
B,  respectively,  using  Seeding  Spectrum  II.  Neutral  seedings 
are  not  known  but,  as  before,  produced  negligible  change. 

The  charged  seedings,  however,  produce  more  improvement  in 
both  fogs  with  both  seeding  concentrations  than  was  previously 
produced  with  the  larger  sized  seeding  drops.  Compared  to 
Figures  8 and  10,  Figures  13  and  14  show  that  more  fog  liquid 
water  is  removed  in  both  fogs.  Practically  all  of  the  fog 
LWC  is  removed  in  Fog  A.  Unfortunately,  a considerable  amount 
of  the  treatment,  which  contains  the  fog  liquid  water,  is  left 
behind.  The  improved  visibilities  created  by  Spectrum  II 
suggests  that,  although  approximately  the  same  total  amount 
of  liquid  water  remains  in  the  air  (compared  to  Spectrum  I), 
water  has  been  transferred  into  larger,  less  visibility- 
restricting  sizes.  A similar  change  occurs  in  Fog  B,  but  with 
even  more  total  liquid  remaining  (com paring  Figure  14  to 
Figure  10).  Obviously,  a relatively  small  concentration  of 
seeding  liquid  can  result  in  the  collection  of  most  of  the 
fog  drops.  However,  unless  the  growth  of  the  seeding  parti- 
cles results  in  final  sizes  large  enough  to  fall  from  the 
fog,  the  seeding  material  itself  will  remain  to  limit 
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Figure  11.  Visibility  as  a function  of  time  (min)  and  height  using  Seeding 
Spectrum  II  in  Fog  A. 
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Figure  12.  Visibility  as  a function  of  time  (min)  and  height  using  Seeding 
Spectrum  II  in  Fog  B. 
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Figure  13.  For  alternate  model  levels,  the  change  in 
fog  and  treatment  LWC  for  experiment  depicted  in 
Figure  lib  (Fog  A with  charged  Seeding  Spectrum  II 
seeding  of  0.2  gm  m-^) . 
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Figure  14.  For  alternate  model  levels,  the  change  in 
fog  and  treatment  LWC  for  experiment  depicted  in 
Figure  12b  (Fog  B with  charged  Seeding  Spectrum  II 
seeding  of  0.2  gm  m~3) . 
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visibility  somewhat.  Unlike  the  earlier  study  of  an  applied 
electric  field  in  which  visibility  improvement  resulted 
primarily  from  removal  of  fog  water,  seeding  with  charged 
drops  produces  visibility  change  both  from  this  removal  and 
from  a movement  of  water  from  the  fog  to  the  treatment 
spectrum.  Because  most  seeding  drops  are  significantly 
larger  in  size  than  the  fog  droplets,  movement  of  water  from 
small  visibility-restricting  sizes  to  larger  less-restricting 
radii  is  significant. 

Subsequent  reductions  in  treatment  drop  size  caused 
maximum  visibilities,  in  both  fogs,  either  to  decrease  or 
remain  about  the  same.  In  general,  although  more  and  more 
fog  water  was  collected,  more  total  water  remained  behind. 

And,  because  this  remaining  water  is  associated  with  average 
sizes  smaller  than  those  associated  with  larger  initial  treat- 
ment drops,  visibility  decreases.  It  would  thus  appear  that, 
with  the  concentrations  that  can  be  produced  by  current 
research  aircraft,  a treatment  distribution  averaging  between 
10  and  15  ym  in  radius  is  ideal. 

c . Experiment  Set  III  - Variation  of  Treatment  Concentration 

It  was  earlier  evident  that  a treatment  concentration 

_3 

of  0.2  gm  m produced  better  results  than  half  that  amount. 
With  the  idea  that  more  is  always  better.  Fogs  A and  B were 
seeded  with  ever  larger  concentrations.  One  would  expect, 
however,  that  a point  would  be  reached  beyond  which  increased 
seedin’g  would  prove  detrimental.  Such  was  shown  to  be  true 
for  hygroscopic  seeding  of  warm  fogs  (see,  e.g..  Tag  et  al . , 
1970  or  Silverman  and  Kunkel , 1970).  Although  the  basic 
mechanism  (condensation)  for  that  technique  is  different, 
the  effect  caused  by  over-seeding  is  the  same  --  too  many 
particles  competing  for  the  available  water,  growing  too 
little,  and  remaining  in  the  fog. 
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Figure  15  summarizes  final  visibility  (after  15  min)  as 
a function  of  initial  treatment  concentration.  Both  Fogs  A 
and  B and  Seeding  Spectra  I and  II  are  considered.  Visibil- 
ity at  the  lower  of  the  two  seeding  levels  (135  m,  level  9) 
is  plotted.  Obviously,  an  optimum  concentration  for  one  fog 
and  seeding  spectrum  is  not  the  same  as  for  others.  In 
general,  as  the  mean  radius  of  the  seeding  spectrum  decreases, 
smaller  concentrations  become  optimum.  It  must  be  remembered, 
however,  that  the  operational  situation  described  earlier 

could  only  produce  average  concentrations  on  the  order  of 
_ 3 

0.1  gm  m . As  a result,  unless  a particular  area  were 
repeatedly  treated  with  multiple  seeding  passes,  the  threat 
of  over-seeding  is  not  a concern.  Even  small  pockets  of 
initially  large  concentrations  left  in  the  seeding  wake 
should  not  remain  for  very  long.  As  is  evident  from  Figure 
15,  for  the  seeding  spectrum  currently  in  use  in  the  field 
(Spectrum  I),  seeding  rates  could  be  increased  up  to  eight 
times  with  ever  better  results. 

d . Experiment  Set  IV  - Variation  of  Seeding  Drop  Charge 

A comparison  of  the  collection  kernels  in  Figures  3b 
and  c,  for  15  and  30%  of  the  Rayleigh  Limit  respectively, 
shows  the  effect  of  increased  charge.  Equation  (2)  indicates 
that  a doubling  of  charge  results  in  a 1.74  increase  in 
collision  efficiency,  and  thus  the  sweepout  rate.  Several 
experiments  were  conducted  to  determine  how  this  increase  is 
transTated  into  visibility  improvement. 

An  initial  charge  of  30%  of  the  Rayleigh  Limit  was  tried 
with  both  seeding  spectra  in  both  fogs.  It  was  found  that 
a doubling  of  charge  results,  generally,  in  a doubling  of 
visibility  improvement.  No  significant  trends  were  discov- 
ered for  either  fog  or  seeding  spectrum.  The  charge  increase 
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Figure  15.  Final  visibility  as  a function  of 
ment  concentration  for  both  seeding  spectra 


initial  treat- 
in  both  fogs. 
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gave  a slightly  larger  percentage  of  improvement  in  the 
lower  fog  levels.  Also,  those  situations  which  had  resulted 
in  large  improvements  for  the  15%  seedings  gave  smaller 
percentage  increases  than  those  situations  which  had  respon- 
ded poorly.  Technical  considerations  in  the  field  will 
surely  dictate  limits  to  the  eventual  seeding  drop  charge, 
but  theory  suggests  that  charge  be  maximized  as  fully  as 
possi bl e . 
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7.  AN  EVALUATION 


From  this  study  alone,  a complete  evaluation  of  the 
feasibility  of  using  charged  water  drops  to  clear  fog  is 
impossible.  The  one-dimensionality  of  the  model  has  precluded 
processes  which  would  tend  to  limit  clearing  potential.  In 
addition,  there  are  currently  no  field  test  data  with  which 
to  compare  the  results.  A comparison  to  the  earlier-mentioned 
study  by  Carroz  ejt  aJL  is  not  feasible  because  of  the  differ- 
ent nature  of  that  study.  In  addition  to  using  monodisperse 
spectra,  hygroscopic  growth  was  included  in  the  growth  of 
drops  charged  completely  to  the  Rayleigh  Limit. 

One  comparison  which  is  of  significance  relates  to 
earlier  numerical  studies  of  hygroscopic  seeding  of  warm  fogs. 
Of  these  studies,  one  was  made  using  basically  the  same  model 
as  that  used  in  this  study  (see,  e.g.,  Tag  e_t  aj_.  ( 1 970  ) or 
Tag  (1971)).  Condensation,  rather  than  collection  growth, 
was  the  important  mechanism  in  the  hygroscopic  simulations. 

To  allow  an  accurate  comparison  between  hygroscopic  and 
charged  drop  seedings,  several  hygroscopic  seeding  experiments 
were  repeated  using  treatment  concentrations  and  spectra 
duplicating  those  used  in  this  study.  As  in  the  earlier 
hygroscopic  seeding  studies  and  that  described  by  Carroz 
et_  aj_.  , a mixture  of  urea,  ammonium  nitrate,  and  water  is  the 
hygroscopic  solution  being  simulated. 

Figures  16  and  17  represent  hygroscopic  seedings  of  Fogs 

- 3 

A and  B using  a concentration  of  0.1  gm  m with  Seeding 
Spectrum  I.  They  can  be  compared  directly  to  Figures  7b  and 
9b,  respectively.  Note  that  the  abscissa  is  now  logarith- 
mically, rather  than  linearly,  scaled.  Obviously,  the  model 
predicts  that  hygroscopic  seedings  will  result  in  much  more 
visiblity  improvement.  VIF's  are  one  to  two  orders  of 
magnitude  greater  than  those  from  the  charged  seedings. 
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HYGROSCOPIC  SEEDING  0.1  gm  m 
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VISIBILITY  (m) 

Figure  16 . Visibility  as  a function  of  time  (min)  and  height  using  a hygroscopic 
Seeding  Spectrum  I in  Fog  A. 
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Figure  17.  Visibility  as  a function  of  time  (min)  and  height  using  a hygroscopic 
Seeding  Spectrum  I in  Fog  B. 


In  contrast  to  the  potential  suggested  by  numerical 
modeling  (see  also  Silverman  and  Kunkel  (1970)),  field  studies 
of  hygroscopic  seeding  have  produced  inconsistent  results. 
Under  relatively  stable  maritime  conditions,  St.  Amand  et  al . 
(1971)  produced  repeatable  clearings  in  seven  out  of  seven 
experiments;  in  later  tests  at  the  same • 1 ocati on , however, 
Hindman  ejt  a_l_.  ( 1 973)  reported  that  three  out  of  six  experi- 
ments raised  visibilities  above  landing  field  minimums.  In 
a separate  set  of  experiments,  Smith  et_  aj_.  ( 1970  ) observed 
clearings  in  three  out  of  ten  tests. 

Thus,  field  studies  of  hygroscopic  seeding  have  not  borne 
out  the  potential  suggested  by  numerical  modeling.  Conse- 
quently, the  above  numerical  comparisons  of  hygroscopic  and 
charged  seedings  should  hold  some  significance.  They  suggest 
that  the  treatment  concentrations  and  charges  that  have  been 
simulated  in  this  paper  are  not  adequate  for  clearing  fog. 
Furthermore,  it  is  concluded  that,  unless  concentrations  and 
charges  can  be  very  greatly  increased,  charged  drop  seeding 
may  not  be  a viable  fog  dissipation  technique. 
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8.  CONCLUSIONS 


The  following  specific  conclusions  were  reached  from  the 
modeling  experiments: 

(1)  Visibility  improvement  is  closely  linked  to  the 
size  of  the  fog  droplets.  For  the  same  treatment  spectrum, 
visibility  improvement  decreases  with  the  average  fog  droplet 
size. 

(2)  For  the  same  amount  of  seeding  water,  a treatment 
spectrum  with  an  average  size  between  10  and  15  ym  is  ideal. 

A distribution  having  too  much  water  concentrated  in  radii 
less  than  10  ym  results  in  drops  which,  after  growth,  have 
fall  velocities  that  are  too  small.  A distribution  with  too 
many  large  drops  --  although  each  drop  is  quite  efficient  as 
a collector  --  does  not  possess  the  number  density  required 
to  optimally  improve  visibility. 

(3)  For  an  applied  electric  field,  it  was  shown  that 
visibility  improvement  results  primarily  from  droplet  fallout 
caused  by  spectral  widening.  With  charged  drop  seeding, 
improvement  results  both  from  a removal  of  foq  water  (to  the 
ground)  and  from  a transfer  of  water  from  the  fog  spectrum 

to  the  larger  treatment  drops. 

(4)  Visibility  improvement  increases  with  seeding  rate. 
For  a given  seeding  spectrum,  a theoretical  point  of  over- 
seeding is  possible  but  unlikely  with  concentrations 
generated  by  current  airborne  systems  capable  of  dispensing 
charged  water  drops. 

(5)  Visibility  improvement  increases  with  seeding  drop 
charge.  This  charge  should  be  maximized  as  fully  as  possible. 
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The  following  conclusions  were  reached  by  comparing 
results  of  this  study  to  results  of  simulations  of  hygroscopic 
seeding: 

(1)  The  charges  and  treatment  concentrations  that  have 
been  simulated  in  this  paper  would  not  be  adequate  for 
clearing  fog. 

(2)  Unless  charges  and  seeding  concentrations  can  be 
very  greatly  increased,  charged  drop  seeding  is  probably  not 
a viable  fog  dissipation  technique. 

It  is  hoped  that  results  from  future  field  experiments 
will  allow  comparisons  which  test  these  latter  two  conclusions. 
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